Excitations in a non-equilibrium Bose-Einstein condensate of exciton-polaritons 
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We have developed a mean-field model to describe the dynamics of a non-equilibrium Bose- 
Einstein condensate of exciton-polaritons in a semiconductor microcavity. The spectrum of ele- 
mentary excitations around the stationary state is analytically studied in difi'erent geometries. A 
diffusive behaviour of the Goldstone mode is found in the spatially homogeneous case and new 
features are predicted for the Josephson effect in a two-well geometry. 
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After a few decades of impressive efforts on a vari- 
ety of different systems such as bulk cuprous oxide Q 
and coupled quantum wells first observations of 
Bose-Einstein condensation (BEC) of excitons in solid 
state systems have been recently reported in a gas 
of exciton-polaritons Q and immediately confirmed by 
other groups 0, . The system under investigation con- 
sists of a semiconductor microcavity containing a few 
quantum wells with an excitonic transition strongly cou- 
pled to the cavity photon mode. In this strong coupling 
regime, the basic excitations of the system are exciton- 
polaritons, i.e. linear superpositions of a quantum well 
exciton and a cavity photon. As compared to other exam- 
ples of BEC, namely in liquid ^He and ultracold atomic 
gases, the main novelty of the present polariton system 
is its intrinsic non-equilibrium nature due to the finite 
lifetime of polaritons. The condensate has in fact to 
be continuously replenished from the relaxation of op- 
tically injected high energy excitations (e.g. free carri- 
ers or hot polaritons), and its steady state results from 
a dynamical equilibrium between pumping and losses. 
This makes the present system a unique candidate for the 
study of the BEC phase transition in a non-equilibrium 
context. Recent theoretical work 6] has suggested that 
the non-equilibrium condition is responsible for dramatic 
changes in the dispersion of low-lying excitations of inco- 
herently pumped polariton condensates: the sound mode 
of equilibrium condensates is replaced by a diffusive mode 
with fiat dispersion, as it typically happens in coherently 
driven pattern forming systems, such as Benard cells in 
heat convection [7] or optical parametric oscillators 

The present Letter is devoted to the development of a 
simple and generic model of a non-equilibrium conden- 
sate which does not involve the microscopic physics of the 
polariton, and can be used to describe the dynamics inde- 
pendently of the details of the specific pumping scheme. 
Our model is inspired by classical treatments of laser op- 
eration 0, and closely resembles the generic model of 
atom lasers developed in In this way, we are able 

not only to confirm the conclusions of Ref. Q but also 
to analytically relate the elementary excitation spectrum 
to experimentally accessible quantities. The same model 



is then applied to the Josephson effect 12, in a 
system of two weakly coupled polaritonic condensates: 
predictions are given for the frequency and the intrinsic 
damping rate of Josephson oscillations, and overdamped 
behavior is anticipated in the case of strong damping. 

The experimental scheme used to create the polariton 
condensate is sketched in FigU^: under a continuous- 
wave high energy illumination, hot free carriers are gen- 
erated in the semiconductor material forming the micro- 
cavity. Their cooling down by phonon emission leads 
the formation of a incoherent gas of bound excitons in 
the quantum wells, which eventually accumulate in the 
so-called bottleneck region above the infiection point of 
the lower exciton-polariton (LP) branch 14]. Polariton- 
polariton collisions are then responsible for the (gener- 
ally slower) scattering of polaritons from the bottleneck 
region to the bottom of the LP branch. For high enough 
polariton density, Bose stimulation of scattering into the 
lower part of the LP can take place 15]. When the stimu- 
lated scattering rate overcomes losses, the polariton field 
becomes coherent, and a Bose-Einstein condensate ap- 
pears [lit. 

Our model is based on a mean-field description of the 
coherent polariton field in terms of a generalized Gross- 
Pitaevskii equation (GPE) for the condensate macro- 
scopic wavefunction including loss and amplification 
terms 
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As we are interested in the lowest part of the LP dis- 
persion, a parabolic dispersion can be taken for the po- 
laritons, with an effective mass m^p. The strength of 
polariton-polariton interactions within the condensate is 
fixed by the coupling constant g, while 7 is the polari- 
ton damping rate at the bottom of the band. Assuming 
that relaxation processes are fast enough to ensure local 
equilibrium in the polariton reservoir in the bottleneck 
region, and that all coherences between the reservoir and 
the condensate decay on a fast time-scale compared to 
the condensate dynamics, the state of the reservoir is 
fully determined by its polariton density np(x). The 
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amplification rate of the condensate due to stimulated 
scattering of polaritons from the reservoir is a monotoni- 
cally growing function R (n/j) of the reservoir density nji. 
Interactions between the condensate and the reservoir po- 
laritons are modelled by the interaction constant g, gen- 
erally different from the condensate one g. With respect 
to single-mode theories [13, [13 , the generalized GPE has 
the important advantage of fully taking into account the 
multi-mode nature of the spatially extended polariton 
field. Moreover, unlike kineti c ap proaches based on the 
Boltzmann equation 
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able to describe the coherent dynamics of the conden- 
sate. These issues are essential in view of a study of the 
elementary excitations of the condensate. 

The evolution equation ^ for the condensate has to be 
coupled to an equation for the density nii{x) of reservoir 
polaritons: 

dnji 



dt 



= P-lBnR-R {nn) {x)\' + DV'n.R. (2) 



into HI) and 
is present V'o 



Pth)/l, and the oscillation frequency 
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of the macroscopic wavefunction is /iy = M + 'ign%! with 
As usual 
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22( 1 ■ the elementary excitations spec- 
trum around the stationary state of the system can be 
obtained by a linearization of the motion equations ([T][2]) 
around the steady state solution ([3]|1]). Thanks to the 
translational invariance of the system, the fluctuations 
can be decomposed in their Fourier components: 



V'(r,t) 
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Introducing the fluctuation vector = (ui^jVi^^w^) 
and substituting this expansion into the motion equa- 
tions, one is immediately led to the eigenvalue equation 
-^^k^k — LoUk defining the elementary excitations, where 
the generalized Bogoliubov matrix £k is 
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Polaritons are pumped in the reservoir with a rate P and 
relax at a rate jh- The spatial hole-burning effect due to 
the scattering of reservoir polaritons into the condensate 
is taken into account by the R{nj^) \tp\'^ term; spatial dif- 
fusion of reservoir polaritons takes place with a diffusion 
constant D. 

The steady state of the system under a continuous- 
wave and uniform pumping P can be obtained by substi- 
tuting the ansatz 

^|;ix,t) - e-^^-*Vo (3) 
nR{x,t) = n%. (4) 

For small values of P, no condensate 
and the reservoir density is a linear 
function of the pump intensity — P/jr. This solu- 
tion is dynamically stable as long as the amplification 
rate is not able to overcome the losses, i.e. R{n^) < 7. 
The threshold P = Pm corresponds to the value v}^ for 
the reservoir density, which guarantees equilibrium be- 



tween amplification and losses R{n*l^) = 7. When the 



pumping rate P is increased above the threshold, the so- 
lution -00 = becomes dynamically unstable and a con- 
densate appears. Stationarity imposes the net gain to 
vanish, which fixes the reservoir density to the equilib- 
rium value rin = n*o . The condensate density grows as 



(7) 

Here a = P/Pth — 1 is the relative deviation from the 
threshold pumping intensity and the dimcnsionless coef- 
ficient (3 — ti^rR' ip^ji) / R {^'r) characterizes the depen- 
dence of the amplification rate on the reservoir density, 
and 77 = l + a/3. The standard Hartree-Fock value g = 2g 
has been taken. Quite remarkably, the excitation spec- 
trum does not depend on the actual value of the scat- 
tering rate R of the reservoir into the condensate mode: 
only the effective exponent (3 and the relative pumping 
rate a do matter. Of course, the threshold value Pth of 
the pumping intensity does depend (in an inversely pro- 
portional way) on R. 

, relaxation from excited free carriers. 
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FIG. 1: Top (a) panel: sketch of the pumping and con- 
densate formation scheme, (b-e) panels: real and imagi- 
nary part of the excitation spectrum of a homogeneous po- 
lariton condensate as a function of the wavevector k (in 
units of the healing length ^ = \/h/ {rriLP m))- Parameters: 
7ii/7 = 5 and a = (5 = 7/^1 = 1 (b,c); 7^/7 = 1 and 
a — 0.5, /3 — ^/fi = 'Jr/i^ = 1 (d,e). Rescaled diffusion 
constant DniLp/h ~ 5 x 10~*; the spectra are indistinguish- 
able from the D — case. Dashed lines in (b,d): Bogoliubov 
dispersion Q at equilibrium. 

Typical examples of the dispersion of elementary ex- 
citations around the Bose condensed state are shown in 
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Fig. [Ijb-d). As the onset of Bose condensation cor- 
responds to a spontaneous breaking of a U{1) symme- 
try, all spectra involve a Goldstone branch whose dis- 
persion uJn{k) tends to in the long-wavelength fc ^ 
limit 0, UM ■ Physically, this mode can be understood as 
a slow rotation of the condensate phase across the sam- 
ple; the generator (1,-1,0)^ of global phase rotations 
is indeed an eigenvector of £k=o with a vanishing eigen- 
value. 

Let us analyze the different cases in more detail, start- 
ing from the physically most relevant one 7^ ^ 7 where 
the state of the reservoir always remains very close to 
its stationary-state for the given density of condensate 
polaritons. The dispersion for this case is shown in 
Figdjbc): in stark contrast with the linear dispersion 
of the propagating sound mode in equilibrium Bose- 
Einstein condensates pH . , the Goldstone mode (in- 
dicated as -I- in the figure) shows here a diffusive and 
non-propagating behavior at low k. The real part is dis- 
persionless and equal to zero, while the imaginary part 
starts from zero in a quadratic way. This result is in 
agreement with the prediction of the recent work Q 
where these issues were addressed starting from a very 
specific microscopic model of the polaritons, and sug- 
gests that the diffusive behaviour of the Goldstone mode 
is indeed a generic fact in non-equilibrium phase transi- 
tions not only under a coherent pumping as in pattern 
forming systems 0, but also in the case of incoher- 
ent pumping. Note that this diffusive behaviour is in 
no way due to the spatial diffusion of reservoir polari- 
tons, and would be present even in the absence of spa- 
tial diffusion. The value D = 5cm^/s actually chosen in 
the figures is inspired by recent experimental studies on 
CdTe quantum well samples f23|, and corresponds to a 
very small value of the dimcnsionlcss diffusion constant 
i) ^DniLp/h w 5 X 10""*. 

An analytical explanation of this behavior is readily 
obtained by adiabatically eliminating the dispersionless 
and strongly damped reservoir mode {R in the figure) 
whose imaginary part is close to (1 -I- a/?) 'Yr. Taking for 
simplicity D — 0, this leads to the following dispersion 
of the two branches of condensate excitations: 
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UJ±{k) = ± \juJBog{kf - V 
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where ujboq is the usual Bogoliubov dispersion of dilute 
Bose gases at equilibrium 
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The non-equilibrium nature of the system is quantified 
by the effective relaxation rate 



whose value tends to when the threshold is approached 
a > and saturates to 7 for large a 1. The first -I- 
branch of ([8|) is the Goldstone branch which corresponds 
for small k values to a slow rotation of the condensate 
phase, while the — branch corresponds to slow modula- 
tions of the condensate density; for low k values, these 
are damped at the finite rate F defined in (fTO]) . From 
([8|) it is immediate to obtain a prediction for the width 
Afc of the region in k space in which the dispersion of 
the Goldstone mode is fiat Re[wG(fc)] — 0. Outside this 
region, i.e. for k ^ Afc, the ± modes recover the stan- 
dard Bogoliubov modes of an equilibrium condensate. 
In order for a sound-like branch to be observable, one 
however needs that Afc ^ 1/^, ^ = {rriLP /^) being 
the healing length of the condensate. Elementary ma- 
nipulations show that this condition is actually met if 
72«5nOj(l + a/3)7fl//3. 

More complex features are observed in the case where 
7i? and 7 have comparable magnitudes [26^ and the reser- 
voir can no longer be adiabatically eliminated [Fig[TJde)]. 
The most significant feature is the dynamical instability 
Im[w] > of the DI branch that is visible in panel (e) for 
fc^ < 1.2: the homogeneous state ([3]|4]) is no longer dy- 
namically stable and a spatial modulation has to appear 
in the steady-state density profile of the condensate. The 
origin of the instability can be traced back in the repul- 
sive interaction between the condensate and the reservoir 
polaritons: because of the mean-field interaction poten- 
tial gnp, a local depletion of the reservoir density npi(x) 
creates a potential well which attracts the condensate 
polaritons making the local reservoir density drop even 
more. This scenario is only possible for sufficiently slow 
reservoir relaxation rates 7^, < 27/ (1 -f a/?), so that the 
effect of repulsive interactions between the reservoir and 
the condensate is able to overcome the tendency of den- 
sity modulations to relax down as previously discussed. 

As a final point of the Letter, we briefly address the 
effect of the non-equilibrium condition on the Joseph- 
son oscillations between a pair of spatially separated po- 
lariton condensates trapped in adjacent wells. Following 
classical work on the Josephson effect [ill . [l2l . 21 [, the 
dynamics can be described by the following set of equa- 
tions 

^^ = -^^3-, + U IV^.f + [R{n,) - 7] V-, (11) 
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T = al3-f/{l + al3), 



(10) 



for the amplitude V'j of the two condensates (j = {1, 2}), 
and the reservoir densities Uj. These equations can be 
derived by projecting the GPE onto the localized wave- 
functions (f>j describing the ground state of each well, and 
assuming these not to be distorted by the interactions. 
The total polariton wavefunction then reads: '0(r) = 
ijji 0i(r)-|-7/)2 02 (r), the charging energy [/ = g Jdr\(f>j\^, 
and the hopping energy is related to the polariton ffux 
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through a surface separating the wells by 2l|, [2 



J 



da (0* dz4>2 - (t>*2 dz4>i)- 



(13) 



The standard normalization /dr|(/)jp = 1 has been as- 
sumed here. Inserting typical values of 2 /zm for the 
well size, a separation between the wells of 1 /xm, and 
a well depth of 3meV, one obtains that J = 0.1 meV 
and U = 0.03 mcV for a polariton nonlinearity g = 
0.015 meVAim^. 
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FIG. 2: Time evolution of the population A'^i after an excita- 
tion sequence of duration 'yTexc ~ 1 (dashed vertical lines). 
Upper (a) panel: Josephson oscillations (J = 7/2). Lower (b) 
panel: overdamped relaxation (J = 7/2O). Other parameters: 
U = J/10, 7fl = 10 7, Ro = 7, Po/Pthr = 2, AP/Po = 0.5. 

Restricting again to the most significant 7_r 3> 7 case, 
the frequency of the small amplitude Josephson oscilla- 
tions around the stationary state with Ni^2 = 1^^1,21^ = N 
atoms per well can be obtained by simply replacing the 
expression ujj — y/UiNU -h J) to the Bogoliubov fre- 
quency LOBogik) in ([8]). Examples of the different regimes 
for the dynamics around the stationary state are shown 
in Figl^l starting from the steady-state under Pi^2 = 
Pq, the pumping intensity in each well is modulated to 
Pi^2 = -Poi AP for a short time interval < t < T^xc and 
then brought back to Pi 2 = -Pq- The system dynamics 
is followed on the mode populations A''! 2- If the effective 
damping rate T is smaller than ujj (upper panel), the 
behaviour closely ressembles the usual Josephson oscil- 
lations as observed in atomic Bose-Einstein condensates 
in [l3| : the only differences are the intrinsic damping rate 
r, and the slight frequency shift predicted by dH). On the 
other hand, if F > (lower panel), Josephson oscilla- 
tions are replaced by an exponential relaxation back to 
the stationary state; the two modes at u}± predicted in ([8]) 
appear in the relaxation dynamics as two well-separated 
exponentials. 

In summary, we have developped in this Letter a 
generic model for the coherent dynamics of a polariton 
Bose-Einstein condensate. The intrinsic non-equilibrium 



nature of the system is taken into account by means of a 
generalized Gross-Pitaevskii equation including loss and 
amplification terms. This model is used to study the el- 
ementary excitations around the stationary state of the 
system. In a spatially homogeneous geometry, a diffu- 
sive Goldstone mode is found, and novel features are pre- 
dicted for the Josephson effect in two- well systems. This 
formalism will be of great utility in the study of the com- 
plex structures which appear in spatially inhomogeneous 
condensates thanks to the interplay of condensation and 
losses (see, e.g., [29|). 
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